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Abstract 

The problem of neutron scattering by the single magnetic atom is theoretically considered in the second 
order perturbation theory. It is demonstrated that elastic scattering of unpolarized neutron by magnetic 
atom is skewed, i.e. contains the term with the symmetry of mixed product of the atom magnetic moment 
and wave vectors of incident and scattered neutrons ([k x k'] ■ h). The problem of dynamical diffraction 
of unpolarized neutrons by the perfect ferromagnetic crystal is investigated. The case is considered when 
the Bragg condition is satisfied for two reciprocal lattice vectors. In this case neutron skew scattering 
manifests itself as the dependence of the diffracted beam intensity on the sign of the crystal magnetization. 
The calculations have been done for the diffraction of unpolarized neutrons by Dy crystal. In case of the 
diffraction by Dy the change of the intensity under the magnetization reversal achieves 50%. 

Introduction 

As it is well known a neutron is the neutral particle and there are no spin-orbit interaction and 
Lorenz force for it. This leads to absence of Hall-like effects for neutrons moving in uniform magnetic 
field. Magnetic field influences on neutron behavior due to existence of a magnetic moment of the 
particle. Magnetic dipole interaction is described by the operator —fi(aB) [1], where fi is the neutron 
magnetic moment, and B is the induction of magnetic field, a is Pauli matrices vector. It is well known 
that the interaction of this kind leads to Hall-like effect in systems with non-uniform magnetic field. 
This phenomenon, so-called topological Hall Effect (THE), has been theoretically predicted [2,3] and 
experimentally observed [4] for conduction electrons in solids. In this paper effect similar to the topological 
Hall Effect is theoretically investigated for neutrons. 

In the case of electrons Pauli term describes the exchange interaction of conduction electrons with 
localized ones. It has the form —J(aM), where J is an exchange constant, M is the unit vector co 
directed with magnetization. This term can be reduced to the effect of fictitious Lorenz force which 
is non zero only in the systems with non-coplanar and at least two-dimensional magnetization spatial 
distributions. It is naturally to expect appearance of the effect for neutrons similar to THE in the same 
magnetic systems. However in contrast to the exchange interaction depending on the magnetization of the 
media the magnetic dipole interaction of neutrons contains the induction of the magnetic field. Therefore 
exactly the magnetic induction should satisfy the criteria mentioned above. The simplest example of three 
dimensional non-coplanar magnetic induction spatial distributions is the field of the point dipole. The 
scattering of unpolarized neutrons by a magnetic dipole in the first order Born approximation does not 
show any Hall-like effects [1]. But beyond the first order of the perturbation theory new effects appear. 
Some kind of Hall Effect for polarized neutrons was firstly described by Maleyev and Toperverg [5,]. It 
has been shown that inelastic scattering cross section of polarized neutrons contains the term (P[k x k ]), 
where P is the polarization of incident neutron beam, k, k' are the wave vectors of the incident and 
scattered beams. This effect was observed experimentally by Okorokov, Gusakov, Otchik and Runov [6]. 
In the present paper it is shown that the elastic scattering cross section of unpolarized neutrons by a 
single magnetic dipole can contain the term (S[k x k']), where S is the dipole magnetic moment. This 
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term also arises in the second order perturbation theory. Such a peculiarity of cross section can be named 
as "neutron skew scattering" by analogy with skew scattering of electrons in ferromagnets [7] . 

In the ferromagnetic crystals magnetic dipoles form the regular lattice. The peculiar properties 
of the neutron scattering by the single dipole should reveal itself in the diffraction by such a crystal. 
In the present paper the dynamical diffraction of unpolarized neutrons by the perfect ferromagnetic 
crystal is investigated. In contrast to previous works on this subject [8-10 and many other] "three wave" 
approximation is considered, i.e. the approximation in which the Bragg condition is nearly satisfied for 
two reciprocal lattice vectors. Only in this case non-coplanar character of the microscopic magnetic field 
plays the role and the neutron "skew" scattering makes the contribution oc (M[k x k']) to the diffraction 
(M is the crystal magnetization). 

The paper is organized as follows. In the second section the neutron cross section by a single magnetic 
atom is calculated in the frame of perturbation theory. In the third section the unpolarized neutron 
diffraction by a perfect ferromagnetic crystal is considered. 

Magnetic elastic scattering of unpolarized neutron by point magnetic dipole 

Consider cold neutron magnetic scattering by an atom. Spin dependent scattering amplitude of such 
a process is given by the formula [11] 

P\k, k') = 2r ol SP(k-k')0 ■ q m ), (1) 

q m =h-(h-e)e, (2) 

e = (k-k')/\k-k'\. (3) 

Here J 1 is a spin matrix, superscript (1) means that the scattering amplitude is calculated in the first 
order of perturbation theory, vq is classical electron radius, 7 = —1.913 is a neutron magnetic moment 
expressed in the units of nuclear magneton, S is the atom spin (spin of the electrons), P(k — k') is a 
magnetic form factor, h is the unit vector along the magnetic moment of atom. In the sake of simplicity 
lets consider the point dipole, i.e. assume P(k — k') = 1. The cross section of unpolarized neutron elastic 
magnetic scattering am\k,k') in the first order Born approximation is given by the well known formula 
[1]: 

a£(k,k') = (r o7 ) 2 S 2 (l-(eh) 2 ) (4) 

The scattering is reciprocal 

a£\k,k')=a£\-k>,-k) (5) 

and does not contain term (h[k x k']). 

In the second order Born approximation the scattering amplitude has the form 

J e - h z q z /2m n + id 
AttH 2 

V m (k-k') = S(a-q m ) (7) 

m n 

e = h 2 k 2 /2m n is neutron energy, m n is neutron mass. This scattering amplitude corresponds to the double 
neutron scattering. Correction to the scattering cross section of unpolarized neutron () is determined by 
the interference of the first and second order of perturbation series for scattering amplitude. 

ag\k,k') = 2TrRe(f 1 \k,kT T PHk,k')). (8) 
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Tr means the trace over spin coordinates, i?emeans real part. It is easy to show that 

Tr(V m (k' - ky T V m {k' - q)V m (q- k)) = i(q ml [q m2 X q m3 }), (9) 

q m i=h-{h-e 1 )e 1 , e l = (k - k')/\k - k'\, (10) 

q m2 = h-(h- e 2 )e 2 , e 2 = [k' - q)/\k' - q\, (11) 

q m3 =h-(h- e 3 )e 3 , e 3 = (q - k)/\q- k\. (12) 

Correspondingly 

a [ l\k, k') = ^i(r o7 ) 3 5 3 / 5(s - tfq 2 /2m n )(q ml [q m2 x q m3 })d 3 q, (13) 

In the coplanar situation, i.e. when (h ■ [k x k'}) — 0, there are no second order correction to the cross 
section dm (k, k') = 0. Indeed the expression (8) in this case transforms to 

a£\k,k')cx([hx(k-k')}-F), (14) 

F = J \{h- e 2 )(h ■ ei) + (h ■ e 2 )(h ■ e 3 ) + (h ■ e 3 )(h ■ e*i)} S(e - h 2 q 2 /2m n )d 3 q. (15) 

Let's k, k' and h lays in the plane (x,y), then only z component of the vector F gives contribution to 
<Jm\k, k'). The z component of F is zero since the products (h-ei) depend only on q\, and do not depend 
on q z . 

For the case of general orientation of k, k' and h the contribution a^m {k, k') to the scattering cross 
section is an odd function of h. Also <Tm\k, k') changes the sign under permutation k — >• k! , k! —t k and 
under permutation k — > —k', k' — > — k. Thus the term <Jm\k, k') has the same properties as the mixed 
product (h ■ [k x k']). 

The dependence of the unpolarized neutron scattering cross section (cr m (fc, k') — am\k,k') + 
<Jm (k, k')) on the angle between wave vectors of incident and scattered neutrons is presented. Vec- 
tors k and k' are perpendicular to the magnetic moment of point dipole. For better representation the 
case is presented of unfeasibly strong magnetic interaction at which dm\k, k') is of order of <Jm\k, k'). 
For real magnetic atoms and wavelengths of cold neutrons the contribution of the second order is very 
small <Jm\k, k')/<Jm\k, k') oc r j/\ ~ 10~ 4 (A is neutron wavelength). 
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(a) (b) 

Fig. 1. Dependence of uploarized neutron scattering cross section (in arbitrary units) on the angle tp 
between wave vectors of incident (k) and scattered (Ar) beams. Two orientations of dipole magnetic 

moment. Polar plot. 



It can be seen from the Fig. 1 (a) the neutrons scatter mainly to the right with respect to incident 
beam and the scattering is indeed skewed. If one reverse the direction if the magnetic dipole the preferred 
direction of the scattering also is reversed with respect to incident vector wave vector. 

Diffraction of unpolarized neutrons by the perfect ferromagnetic crystal 

It is reasonable that the phenomena of neutron skew scattering can be enhanced by means of the 
interference effects. Perfect ferromagnetic crystal represents regular lattice of codirected magnetic dipoles. 
Let's consider dynamical diffraction of unpolarized neutrons by such a crystal. 

At first analyze the results of the previous section to understand the necessary condition for "skew" 
diffraction of neutrons. As it was shown the "skew" scattering appears due to the interference of the once 
and twice scattered waves (see Fig. 2). 
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Fig. 2. Scattering from the state k to the state k! . Interference of the once and twice scattered waves. 

In the case of the scattering by the single dipole all intermediate states q on the energy surface con- 
tribute to the scattering. In the case of diffraction only the certain direction of scattering are permitted. 
Therefore possible intermediate states determined by the Bragg conditions. If the Bragg condition is 
fulfilled only for one reciprocal vector (here this situation named as "two wave approximation" ) it intu- 
itively seems that the double scattering processes through the intermediate state (depicted by the dash 
line on the figure) are not taken into account. Indeed as it has been shown earlier [8] the diffraction in 
this case is reciprocal 

I(k,k',M) = I(k',k,-M), (16) 

I(-k',-k,M) = I{k,k',M), (17) 

and therefore does not contain terms of the from (M ■ [k x k']). In the formulas (16), (17) / denotes the 
intensity of neutrons diffracted in the direction k' {k' — k = g, g is the reciprocal lattice vector). 

Double scattering through the intermediate state appears if one considers the "three wave approxi- 
mation". In this case the Bragg condition is nearly satisfied for two reciprocal vectors (see Fig. 3). 

From the formula (8)it is seen that the contribution to the "skew" scattering is made by the intermedi- 
ate states with vectors q having non-zero projection on the direction of the atom magnetic moment. Thus 
if all the reciprocal lattice vectors participating in the diffraction process perpendicular to the magnetic 
moment of the crystal no skew scattering should appear. 

Let's now investigate the problem of "three wave diffraction" for unpolarized neutrons. In accordance 
with Bloch theorem the wave function of neutron inside the ferromagnet is written as follows 

* £ = J2 ^%e l ^+^ (18) 

i 

Here k is quasimomentum of a neutron in media, gi are the reciprocal lattice vectors, are the two 
component amplitudes. The system of equations for ty gi has the form 

(1 - k 2 / X 2 )^% - J2(N 9 k + M gk )* k gi _ gk = (19) 
k 
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Neutron energy e is related withx by the formulae = ?i 2 x 2 /2m n . N g i are the coefficients of the Fourier 
series of the nuclear potential. 

N gk = ^^ [ d\V nuc {r)eMigk?). (20) 

11 X Junitcell 

N c denotes the number of unit cells per unit of the crystal volume, V nuc (f) is the nuclear potential. 
In the case of slow neutrons Fermi pseudo potential approximation can be used for the description of 
neutron-nuclear interaction. 

Vnuc if) = — V bS(f- fi) (21) 

m„ 

Summation in the formula is made over all atoms in a crystal, b is the coherent nuclear scattering length. 
Here the dependence of the potential on the nuclear spin is not taken into account since the temperature 
is assumed to be not very low. 
Thus 

N gk = bl eMWkRl), (22) 

x i 

where summation is carried over the atoms in one cell, Ri is the position of atoms. 
M g k are the coefficients of the Fourier series of magnetic dipole interaction. 

M gk = ^-^-(a ■ q gk ) ^2 Pi cxpiigkRi) (23) 
x i 

q gk =ti - (h- e gk )e gk , e gk = g k /\g k \ (24) 

pi denotes magnetic form factor of the l-th atom. In the "three wave approximation" vectors g\ and <?2 
meet the Bragg condition and correspondingly amplitudes^Q , and ^ g2 are non-zero. One should 
solve the following system of equations to find them 

(1 - k 2 / x 2 -No - ao0 ■ 9o))*o - (N gl + ax{5 ■ q g i))^ g i - {N g2 + a 2 (<? • q g2 ))^ g2 = 0, 
-(N gl + a o ■ q g i))y g a + (1 - (k + gif/x 2 - N - ai ■ <f ))* 3 i - (N g3 + a 3 ■ <7 33 ))* 9 2 = 
-{N g2 + a 2 ■ q 9 2))^ g o - {N g3 + a 3 ■ g s3 ))* s i + (1 - (fe + 9i? Ix 2 - N - a {S ■ <? ))* S 2 = 

(25) 

Here the spin dependent part of neutron-atom interaction is presented in the detailed form. This equations 
contain four magnetic vectors <fo-3 (<73 = h— (h-e g 3)e g 3, e g3 = (gi — .92) / |.9i ~9i\)- F° r simplicity propose 
that gi is perpendicular to magnetization gi-Lh, then q = qi = h, and there are three different magnetic 
vectors. 

Let's consider the symmetry of the equations (25) with respect to magnetization reversal, i.e. to the 
replacement h — > —h. Denote this operation by T. Terms connected with magnetic interaction changes 
their sign under the reversal . 

fq gi = -q gi . (26) 

In the case when all three vectors lay in one plane the action of the operator T can be compensated 
by the spin coordinate system rotation by ir around the axis perpendicular to the vectors <fo-3- It is 
clear that spin coordinates rotation does not change the intensity of the diffracted beams. Therefore the 
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diffraction of unpolarized neutrons by the perfect ferromagnetic crystal in this case is an even function 
of magnetization and does not contain the contribution of "skew" scattering. 

If <fi_ 3 is non-coplanar the action of the operator T can not be compensated by any spin coordinates 
rotation. Thus the intensity of the diffracted wave can depend on the sing of the magnetization. The 
condition of non-coplanar alignment of the vectors 9*1-3 in case of gi-Lh can be formulated in the form 

(5 2 fe) 2 (fe-[.9ix.g 2 ])^0. (27) 

On the base of the equations (25) the diffraction of unpolarized neutrons by Dy crystal has been 
calculated. The following geometry is chosen. Magnetization of the crystal is parallel the interface (Fig. 
3). The unpolarized neutron beam falls onto the crystal surface with the glancing angle a. The wave 
vector of the incident neutrons ki nc is perpendicular to the crystal magnetization. The vector g\ is 
perpendicular to the crystal surface and also to the magnetization. So one of the diffracted wave comes 
out from the front surface with the wave vector kdiff- The intensity of this diffracted wave is calculated 
for two direction of the magnetization M and —M. The second reciprocal lattice vector §2 is tilted with 
respect to magnetization by the certain angle determined by the Bragg condition. 




Fig. 3. "Three wage" geometry. Bragg case, k is one of the possible neutron wave vectors inside the 

crystal. 

Analytical solution the system (25) is hardly possible since there are 12 neutrons waves with the same 
energy in the crystal, and the dispersion relation arising from equations (25) is the 12-th order equation 
with complicate coefficients. Roots of the dispersion relation have been found numerically. Using the 
roots and corresponding border conditions [8] (the case of semi infinite crystal has been considered) the 
intensity of the diffracted waves has been found. Since the glancing angle is not very small the intensity 
of the reflected wave several order smaller than the diffracted one. 

Dy has the simple hexagonal crystal structure. Vector g\ is codirected with c-axis, and 52 is along 
a-axis. The following parameters of the crystal has been used: magnetic moment per atom fi = 10.8 
li nuc , the coherent scattering length b = —1.4 • 10~ 12 cm, lattice parameters a = 0.35 nm, c = 0.45 nm. 
The wave length is A = 0.25 nm. 



7 



The dependence of the intensity of the diffracted beam on the glancing angle for two magnetization 
direction is presented on Fig. 3. 
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Fig. 4. The dependence of the intensity of the diffracted beam on the glancing angle. The intensity is 

normalized by the intensity of the incident beam. 

As it is seen from the figure the intensity of the diffracted beam depends on the magnetization sign. 
The relative difference of the intensities for two magnetization directions achieves 50% in the narrow 
glancing angle range. Calculations show that the intensity also changes if one interchanges source and 
detector (k inc —> ~kdiff,kdiff — > ~k inc ) or if one send the incident beam from the back side of the 
crystal (the replacement ki nc — > kdiff,kdiff — > ki nc ). Thus one can said the intensity depends on the 
sign of the mixed product ([fc» nc x kdiff] ' h). Such dependence is the signature of the "skew" scattering 
of neutrons by a single magnetic atom. 

In the considered above geometry (Fig. 3) the incident wave vector is laid in the plane perpendicular 
to the magnetic field. If the vector comes out of the plane the Bragg condition for the reciprocal lattice 
vector g 2 can be violated. This leads to disappearing of the effect. Therefore the dependence of the 
diffracted intensity on the sign of the magnetization arises only in the small solid angle region for the 
directions of incident beam. 

Conclusion 

Thus in the paper it is demonstrated that elastic scattering of unpolarized neutron by magnetic atom 
is skewed, i.e. contains the term with the symmetry of ([ft x ft'] • h). Such a term in the scattering cross 
section appears in the second order perturbation theory. In the case of the single atom this contribution 
to cross section is small and is about 10~ 4 . Skew scattering of neutrons manifests itself in dynamical 
diffraction of unpolarized neutrons by perfect ferromagnetic crystal. When the Bragg condition is nearly 
satisfied for two reciprocal lattice vectors the intensity of the diffracted wave depends on the sign of the 
magnetization. The calculations have been done for the diffraction of unpolarized neutrons by Dy crystal. 
The change of the intensity under the magnetization reversal achieves 50% in the small solid angle region 
for the directions of incident wave vector. 
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